Cooling of a particle by coupling to its own reflection 
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We identify and explore a cooling force resulting from the retarded dipole interaction between 
an illuminated particle and its reflection. For a one-dimensional example, we find cooling times of 
milliseconds and limiting temperatures in the millikelvin range. The force, which may be considered 
the prototype for cavity-mediated cooling, may also be enhanced by plasmon and cavity resonances 
at the mirror. 
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The revolution in atomic physics brought about by 
the Doppler coohng of atoms [ij and ions the 
magneto-optical trap Q , and sub-Doppler cooling mech- 
anisms [J, lila, has pr omp ted the proposal of further 
schemes 13,11, B E El E E, Q which might extend 
the ultracold domain to a broader range of species. In 
contrast to the initial reliance upon the scattering force 
of resonant radiation, these more recent schemes have 
largely been based upon the optical dipole force 1^, l3| 
that results when there is a spatially-varying interaction 
with off-resonant illumination. Greatly developed as a 
powerful and flexible means of tweezing atoms and mi- 
croscopic particles the dipole force is at heart con- 
servative and incapable of cooling; when coupled to a 
dissipative component or otherwise invested with a non- 
Markovian character, however, it too can form the basis 
of a cooling mechanism. 

An interesting category of forces to have received much 
recent interest is that in which an atom or ensemble is 
both trapped by, but can also affect, the field within 
an optical cavity, resulting in a cavity-mediated cooling 
force [13, El, m, [13, El- As with the cooling of a parti- 
cle as it moves through a polarization gradient [^ , inten- 
sity variation [3| , speckle field [3] or bichromatic standing 
wave [1] , it is a retarded dependence of the atom-field in- 
teraction upon the atom's position which gives a dissipa- 
tive, velocity-dependent component to the otherwise con- 
servative trapping force. Whereas the retardation in the 
latter cases results from non-adiabatic decay of atomic 
state populations, however, in cavity-mediated cooling it 
is caused by decay of the optical field to which the atom 
is coupled. 

Cavity-mediated cooling may be viewed as a combi- 
nation of two effects: the retarded perturbation by a 
trapped particle of the optical field in which it is trapped, 
as described above; and the amplification of this effect by 
a resonant optical cavity [l8l |. In this paper, we consider 
the simplest case of the former, whereby the perturbed 
field is returned to the atom by a single mirror. This op- 
tical binding phenomenon may be interpreted as 
the classical electrostatic interaction between the dipoles 
induced in the particle and its reflection; or as tweezing of 
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FIG. 1: Retarded binding of a normally-illuminated particle, 
moving with speed v, to its own reflection, depicted (a) in the 
laboratory frame, in which the image lags behind; (b) in the 
rest frame of particle, whereby the 'wake' trails behind. 



the particle by light which the particle itself has focused, 
as shown in Fig. [TJ The finite time taken for light from 
the particle to return via the mirror causes the reflected 
image to trail behind the moving particle, thus provid- 
ing a component of the binding force in the direction of 
the particle velocity. When the sign of this component is 
such as to oppose the particle motion, cooling ensues. 

We analyse the binding between a moving particle and 
its reflection using a similar approach to that previously 
applied to the unretarded binding of two point-like par- 
ticles [2l| . We assume an unperturbed illuminating field 
S'^irA) of wavelength A, particle polarizability a, and 
characterize the round-trip from the particle to its re- 
flection via the mirror by a transit time r and electric 
field propagator (. The electric field f (r^) experienced 
by the particle at position r^(<) and time t is then, to 
lowest order in the velocity v and r and for (a, t|v| <C 1, 
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which corresponds to the result of [2l| with the two par- 
ticles degenerate, modified by the appearance of an ad- 
ditional, velocity-dependent term. With some care over 
the constraints during differentiation, the dipole force ex- 
perienced by the particle may be obtained as in 2l|. A 
series expansion in terms of reveals the leading terms 
for a stationary particle to be the dipole force in the un- 
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perturbed field, then the interaction between two dipoles 
induced by the unperturbed field and, to the same order, 
the dipole force upon the polarized particle due to the 
field propagated from the induced polarization. 

For a one-dimensional geometry, the incident illumina- 
tion combines with its reflection to give the electric field 
a spatial dependence S'^ = £osin2fcx, where k = 27r/A; 
T = 2x/c, and C = —iexp2ikx for propagation from a 
sheet of dipoles a distance x in front of the mirror. With 
a real susceptibility a, the force = ^RelaSS/S*) upon 
the moving particle is therefore 

sin2fcx + 2a (l \ sin^fcx(4 cos^ kx — 1) 



T = 



^iSn k 



2akTv sin4fca 
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The force thus comprises three terms. The first two are 
the dipole force exerted upon the particle by the unper- 
turbed field, and a Dopplcr-shifted optical binding force 
between the particle and its reflection; the Doppler shift 
is manifest here not through the change in frequency but 
through the accompanying change in wavcnumber of the 
components whose superposition yields the field gradi- 
ent. The third term, which depends upon the particle 
velocity, the electric field propagator and the round-trip 
retardation time, is the cooling force, which dominates 
the velocity-dependent part of the second term when the 
distance from the mirror is many wavelengths. The sec- 
ond term is supplemented by a further term in (u/c) due 
to the missing Lorentz component recently noted in [2^ . 

The classical approach is highly instructive, for the 
arbitrary form of r and ( allows the simple geometry of 
Fig.[l]to be extended to more complex arrangements with 
curved mirrors and resonant structures, as well as provid- 
ing a link to the binding and damping of extended or 
highly polarizable particles and micromechanical struc- 
tures [24| . In order to examine the dynamics of a particle 
subjected to the cooling force, however, and to determine 
the limiting temperature at which the cooling is balanced 
by heating from off-resonant scattering, we here outline 
a semi-classical approach [25 1 . 

We examine a one-dimensional geometry in which the 
reflection of a laser-driven atom is imaged back onto it- 
self in a configuration similar to that investigated experi- 
mentally in [26[ , but we consider only a single transverse 
electromagnetic mode corresponding, for example, to the 
use of an optical fiber delay line instead of a free-space 
arrangement. The model thus consists of a single two- 
level atom coupled to the modes of the electromagnetic 
field in the right half space, according to the Hamiltonian 

p2 f 

H ~ hujncr^cr^ + h / Huja' (uj)a(iLj)dijj 
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[x^^ [a+a{uj) - a^uj)a+] duj. (3) 

The first line represents the internal atomic energy given 
by the transition frequency uja, the kinetic energy, and 
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FIG. 2: (a) Spatial dependence of the lowest order friction 
7/m (solid line) for s — 0.1, Ga/{'JmP') = 0.1, a; = 3 m. The 
dashed line indicates the intensity of the pump wave (arbi- 
trary units), (b) Stationary temperature in a tightly confined 
harmonic trap versus trap position for a; = 1 m (dotted), 
X = 3 m (solid), and a; = 10 m (dashed). Here x' is the 
position relative to the nearest field node. 



the mode energies, respectively, and the second line is 
the interaction energy between the atomic dipole and 
the electromagnetic modes with a coupling coefficient g 
which is assumed constant over the range of relevant fre- 
quencies. The density operator p of the atom-field system 
follows the master equation 



dp i , , 



Cp 



(4) 



where Cp is a standard Liouville term corresponding to 
spontaneous atomic decay into free space modes with rate 
2r. For the following discussions we simplify Eqs. ([3]) and 
(HI) by adiabatic elimination of the atomic excited state. 

As a first step, we derive an approximate analytical 
solution for the friction force resulting from the Hamil- 
tonian ([3]). To this end we treat the atomic momen- 
tum p and position x as classical variables and derive the 
Heisenberg equations of motion for the mode operators 
a{uj). These are then approximated by complex numbers 
with initial distribution aito) = A5{lj — ljq) and the sta- 
tionary solution is found to lowest order in the atom-field 



3 



coupling g. We thus obtain the friction force 



with the friction coefRcient 
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7r?ifcV|^p-^sin(4fcx) 



(5) 



(6) 



where A ~ luq — LUa is the pump detuning, r — 2x/c 
is the propagation time of hght from the atom to the 
mirror and back, and we have assumed |A| ^ F. Using 
2Trg'^ ^ ATaa/iirw'^) [13], where aa = 3AV(27r) is the 
atomic scattering cross section and w is the mode waist, 
and introducing the atomic saturation s — (7^|Ap/A^, 
we can rewrite 7 as 
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(7) 



Figure lUa) shows the variation of 7 with atomic position 
in the standing pump wave: regions of coohng (7 > 0) 
alternate with regions of heating (7 < 0). In order to ob- 
tain cooling of a particle by this mirror-induced friction 
force it will therefore be necessary to trap the particle, 
e.g. by a far-off resonant dipole trap, to within A/8; be- 
cause the mean kinetic energy of the trapped particle will 
then be coupled to an equal but uncooled mean potential 
energy, the cooling rate will be halved. For the chosen 
parameters, the theory predicts cooling times (7/m)^^ 
of the order of 3 ms at positions of maximum friction. 

To estimate the steady-state temperature achievable 
using this cooling scheme, we need an estimate of the 
momentum diffusion coefficient. The leading order term 
in the weak coupling limit assumed above is simply given 
by the well-known result for diffusion in the standing 
wave formed by the pump field [2^, i.e., 



D = n^k^Ts 



COS (kx) -\ — sin {kx) 
5 



(8) 



where the factor 2/5 accounts for the spatial distribu- 
tion of circularly polarized spontaneous emission. The 
stationary temperature is hence given by 



ksT ^ - = - 



D h TTw'^ cos'^{kx) 
7 

TTW 



\ sin^(fca;) 



T 2aa 

,2 



sin(4fca;) 



hT- 



(9) 



The final expression of Eq. ^ holds at the position of 
maximum friction. We see that the expected stationary 
temperature is given to within a numerical factor by the 
Doppler limit hV, scaled by the ratio of the cross-section 
TTW^ of the pump beam to that of the atom and with the 
atom-mirror round-trip time r replacing the atomic co- 
herence time 1/F, but is independent of detuning, pump 
intensity and atomic saturation. Figure [D^b) shows the 
stationary temperature of an atom in a tightly confining 
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FIG. 3: (a) Change of temperature versus time for initial 
temperatures of 0.3 mK (solid curves) and 3.1 mK (dashed 
curves). Simulation results are averaged over 10* Monte-Carlo 
trajectories and fast oscillations with the trap frequency have 
been filtered out, straight lines are linear fits. Parameters 
are: mode diameter 1.4 /^m, atomic saturation 0.073, trap 
frequency t'trap = 1-5 MHz, and t = 26.5 ns. (b) Rate of 
temperature change versus initial temperature. Data points 
are simulation results, lines are linear fits. Crosses and solid 
line are for !^trap = 1-5 MHz, and r = 26.5 ns, circles and 
dashed line are for t'trap = 750 kHz, and r = 53 ns. 



harmonic trap versus trap position. For atomic rubidium 
with the same parameters as in Fig.[2)Ja), we obtain tem- 
peratures of the order of 1 mK. While this is higher than 
the Doppler temperature of 140 /iK, the independence 
of detuning means that, for sufficiently off- resonant op- 
eration (|A| ^ lOF), mirror-mediated cooling will be the 
dominant cooling effect. 

We have confirmed the predictions of this simplified 
analytical model by semiclassical Monte-Carlo simula- 
tions of a rubidium atom coupled to a discretized set 
of modes [29| . This model accurately describes quantum 
fluctuations of photon numbers as well as of atomic mo- 
mentum; it is valid for arbitrary coupling strength and 
arbitrary atomic velocity, and an external harmonic trap 
is easily incorporated. However, the numerical model 
does not allow a single trajectory to be followed towards 
a stationary state because of the relatively long cooling 
times (of the order of ms as outlined above) . Instead, we 
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performed the simulations for atoms with various initial 
kinetic energies and monitored the initial change of en- 
ergy with time. Two examples are shown in Fig. ^a). 
For an initially cold sample at 0.3 niK, diffusion dom- 
inates over friction and the temperature increases with 
time. For a larger initial temperature of 3.1 niK, on the 
other hand, friction dominates and a net cooling effect 
is observed; and a sample with exactly the stationary 
temperature remains unchanged over time. Figure El^b) 
summarizes the results for dT/dt obtained from simula- 
tions with varying initial temperatures for two different 
sets of parameters. A linear fit to the initial temperature, 
consistent with a simple random walk model, predicts a 
steady-state temperature of T = 0.61 ± 0.17 mK for trap 
frequency t'trap — 1-5 MHz and delay time t = 26.5 ns, 
and T = 0.31 ± 0.10 mK for ftrap = 750 kHz and r = 53 
ns, where the errors result from the uncertainty of 10"* 
Monte-Carlo trajectories for each set of parameters. For 
trapped particles, our analytical expression ([9]) must be 
modified to account for the period of order t after each 
turning point in the particle's trajectory, during which 
the newly-rcverscd motion is accelerated rather than di- 
minished by the retarded cooling force. For the param- 
eters of Fig. [HI this introduces a further factor of 0.64, 
yielding temperatures of 0.58 mK and 0.30 mK, consis- 
tent with our numerical results. There will be a further 
small correction because the trajectories of the trapped 
particles extend beyond the region of maximal cooling. 

The geometry considered here resembles that in the 
beautiful experiments of [i^, and shares with cavity- 
mediated cooling schemes the requirement for a tightly 
focused image and hence achieves cooling only within a 
small volume. Mirror-mediated cooling, however, per- 
mits a number of possible enhancements. Firstly, the 
mirror may be replaced by an external resonator, such as 
an etalon or multilayer filter, which allows the same retar- 
dation T to be achieved with a greatly reduced mirror dis- 
tance. Material or spatial resonances such as the plasmon 
modes of a micro-structured mirror or antenna might 
achieve this while also enhancing the strength of the scat- 
tered field. As possible examples, we envisage micro- 
scopic antennae [30l or hemispherical, plasmon-resonant 
resonators 3]J, [S^] . These can be produced as arrays so 
that, like a cobbled street for a bicycle, the dissipative 
effect is reproduced over an extended area. 

Two related phenomena are associated with mirror- 
mediated cooling geometries. Firstly, the Casimir-Polder 
force is the effect of retardation upon the van der Waals 
interaction of a vacuum-polarized particle with its reflec- 
tion (33I I : mirror- mediated cooling results instead from a 
small component of the much larger laser-induced force. 
Resistivity in the mirror meanwhile provides an alterna- 
tive dissipative component that results in friction 34 1. 

Although we have here analyzed the longitudinal cool- 
ing of a 1-D geometry, we find that, in three dimensions, 
transverse cooling is also possible. Higher order terms 



can be significant, particularly for more macroscopic par- 
ticles; and we expect that, as with cavity-mediated cool- 
ing, collective effects will be significant for ensembles. 
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